Geometry-induced smoothing of van Hove singularities in capped carbon nanotubes 
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The electronic states of capped semi-infinite nanotubes are studied within the phenomenological 
gauge field-theory model. A single manifold for the description of both the nanotube and the cap 
region (considered as nearly a half of either Ih or I fullerene) is suggested. The wavefunctions and the 

' density of states (DoS) are numerically calculated for both metallic and semiconducting nanotubes. 

The smoothing of van Hove singularities is found and proven analytically. The comparison with the 

• experimental observations is discussed. 
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The electronic properties of capped carbon nanotubes are essential to proposed device applications. The localized 
i— i, eigenstates in a cap region were predicted theoretically in [1]. Experimentally, the signature of resonance states 
localized at the ending caps of a nanotube was observed in the spectroscopic measurements on carbon nanotube 
caps [2] . Another spectroscopic analysis of carbon caps performed in [3] shows spatially resolved electronic properties 
,_h . along capped nanotubes and highlights the van Hove singularities (VHS) in nanotube bulk and resonance states at the 
end. The scanning tunneling spectroscopy (STS) gives information about the local density of states (LDoS). Several 
sharp peaks experimentally observed in [3] were attributed to the VHS. The tight-binding calculations show good 
agreement with the experimental data (see, e.g., [3,4]). 

The aim of our paper is twofold. First, we suggest an appropriate geometry for the shape of capped nanotubes and 
develop a continuum field-theory model to describe the electronic states. Second, we perform numerical calculations 
I \ to obtain the structure of DoS. This study gives us an important information about the influence of the cap region 
on the electronic properties of the nanotube and illuminates the behavior of the VHS. We will consider a closed 
semi-infinite nanotube with a spherical cap. It is suggested that the nanotube region is defect free (contains only 
hexagons) and infinitely long, so that edge effects are negligible. The cap is a half of the spherical (I or Ih) fullerene 
and contains exactly six pentagons. A continuum field-theory model for the description of low-energy electronic states 
in carbon nanostructures of arbitrary geometry was formulated in [5] in the form 
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ON , Here are the conventional Pauli matrices (fi = 1,2), the energy E is accounted from the Fermi energy, the Fermi 
velocity Vp is taken to be one, and the two-component wave function ip represents two graphite sublattices (A and 
B). By using the index k in (1) we take into account electronic states at two independent Fermi wave vectors in 
carbon lattice (so-called "K-spin up" (K') and "K-spin down" (K*) states). Geometry is involved via zweibeins 
e M> = + Qfj, where the spin connection term = (l/8)ui a ^[a a , ap], and (w' 1 )"' 5 are the spin connection 
coefficients: 
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, are the metric connection coefficients. 

In order to take into account pentagons in the cap region, two compensating gauge fields a M and W M are introduced 
in (1). The Abelian field is responsible for the elastic flux due to pentagonal defect [6,7]. Its circulation around a 
single disclination is found to be exactly 27r/6 (the appearance of a pentagon in the honeycomb lattice is equivalent 
to the creation of 60° disclination). The non- Abelian field a allows us to take into account the exchange between 
A and B sublattices in the presence of a pentagon [8,9,5]. Namely, this exchange can be described by using an 
appropriate boundary condition for the K spin part of the four-component spinor wavefunction ip — (ip^ip^) T in 
the form ip(ip + 2ir) = —Tip(<p). Here the holonomy operator T is composed as a product of two operators in the 
form exp(i$Tj), i—2,3 where the isospin Pauli matrices r act on the K part of the spinor components. The general 
consideration takes proper account of the relative placement of pentagons. For even number of defects one has (see [9] ) 

a^=±(iV^+M^). (4) 
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Here the sign plus (minus) is taken for k = (k = K^), respectively, N is a number of defects and M (M = n—m (mod 
3)) depends on the arrangement of pentagons, n and m are the numbers of steps in positive and negative directions, 
respectively. The directions rotated by 2tt/3 are considered to be identical. For (I) fullerenes M = M 2 = ±1 or 
for any two defects, M 4 = for any four defects, and M 6 = —M 2 for six defects. In the case of (Ih) fullerenes M is 
always equal to zero due to the mirror symmetry of the lattice (see [5]). 

In the isospin space, the action of the operator T corresponds to a rotation by 2-kN/A around the second axis 
and then by 2nM/3 around the third axis. The angle between the initial and final vectors gives exactly the phase 
(4). Notice that for odd number of defects the first transformation places the vector on the third axis, so that the 
second transformation does not change the angle (phase) . To simplify a problem, let us introduce an effective field for 
uniformly distributed defects instead of six point-like fluxes due to pentagons in such a way that the phase factors in 
(4) remains the same. Let T$ be the circle on the (hemi)sphere (with the unit radius, for simplicity) which encircles 
the area S including the pole (see fig. 1). Evi< 




FIG. 1. The modelled geometry of closed tube. A contour Fs encircles the surface S. The z-axis is directed downside. 



It is easy to find that the operator of rotation around the second axis in the isospace takes the form exp(ir 2 S'3/2). 
Indeed, we have six defects (each rotating on 27r/4) continuously distributed on the area 27r. The rotation angle is 
$Af = 35/2. The operator of rotation around the third axis can be approximated as M(S) = — M 2 cos(3(S — 27r)/4). 
This agrees with the properties of M described above: M(27r/3) = M 2 , M(47r/3) = 0, M(27r) = -M 2 . The rotation 
angle is written as 



$ T = -^M 2 cos(3(S-27r)/4). (5) 



The total operator T in the form 



T = exp(ir 3 $T) exp(ir 2 $M) (6) 
corresponds to the phase rotation by angle $ which is determined from the equation 

cosfcfS') = cos & M {S) cos $ r (S). (7) 
As a final step, the circulation of the field a along T$ is 

<f a M dx^ =±$(S). (8) 
Jr s 

On the equator of the sphere the circulation is $(27r) = 3tt + 27rM 2 /3. Since the equator is connected with the tube 
region, any contour around the tube encircles the same quantity of the field sources as T 27r . Therefore, finally in the 
spherical region and on the tube we obtain the following conditions: 

/ a v d<p = ±$(S), (9) 

I a v d(p = ±(3tt - 2ttM 2 /3). (10) 
Jo 

Notice that the factor 3ir in (10) can be excluded by the redefinition of the momentum j (see below), so that the 
circulation in the tube region (10) is determined by M 2 in agreement with [10]. For both (Ih) fullerenes and (I) 
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fullcrcncs with Mi = one has $ = <&m = 3.5/2, which agrees with monopole approximation [11,5]. Therefore, one 
can conclude that metallic tubes can be capped by either (Ih) fullerenes (armchair) or (I) fullcrcncs with M2 = 
while the semiconducting tubes are capped only by (I) fullerenes with Mi 7^ 0. 

In the spherical region, the circulation of the field W v has the standard form similar to [5] 



W^dx" 



whereas in the tube region its circulation is equal to — 2tt and may be ignored. 

The continuum description in (1) implies a single manifold with the curvature being a continuous smooth function. 
Therefore we have to propose appropriate geometry for the capped tube: a single manifold £ which reproduces the 
sphere in the cap region and the tube outside. Let us use the following manifold: 



R(p(z) cosy, p(z) sin ip, z), p{z) = Rt\/l — exp(— 2A), (11) 

with 

y I P 

A = f -, a = Rt/Rf, z > -R f , < (p < 2tt, 

R f 

where Rf characterizes the capped area and R t is the tube radius. It is taken into account that Rf can generally be 
different from R t . To simplify the problem, one can approximate S — 2-kK as if the surface would be the hemisphere 
with the radius Rf. The metric tensor is written as 
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The nonzero metric connection coefficients are found to be 

-1 e 2A 
(2(g zz -l) + - 
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+ 1, g vv = p 2 (z), g zv = 0. (12) 



F zz = ~rT~Z (2(<7zz - 1) + — (<7zz I)'), 
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The zweibeins are e\ = ^g zz , e 2 = p(z), and the spin connection coefficients are uj^ = = aR t e~ 2A /(p(z)- s /g zz ), 

so that (1) include the spin connection term 
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Q v = l ^^= . (14) 



Notice that for — Rf < z < the fields are (9), and for z > (10). After the substitution tp = ^ " ) e w the Dirac 
equation (1) takes the form 



aR t e- 2A 
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~ i{ 7k ~ + Wu^ - w * ~ a * ))u = Ev - (15) 

Notice that the momentum j takes half-integer values and enters (15) only by the combination j — a v . We expect 
that in the tube region our model should be identical to that proposed in [10]. Indeed, adding ±3/2 to both j and 
a v , one obtains in (10) f Q n a v d(p = ±27rM 2 /3. Exactly such circulation was suggested in [10] for the effective vector 
potential. 

For large z one has W v = 0, a v = M2/3, y/g zz — 1, p(z) — > Rt, VL V — > 0. In this case, (15) takes a similar to [10] 

form —icrid z ip + (t 2 toV = Eijj, where E = ±\/m 2 + k 2 and the "mass" term m = —(j — A'h/ty/Rt- In this case, the 
general solution to (15) is written as 

^ = c U-L)IE) e ^ (16) 
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with C being a normalization constant. This plane wave solution has nonzero current j z = ipaitp which is valid for 
the infinite tube. In our case, however, the tube is capped and the quasi-particles are not able to leave the tube. For 
this reason, the current has to be zero everywhere. Taking the condition j z — into account we obtain from (16) the 
standing wave solution (ip = ipk + 4>-k) with zero current 

/ _or( cos(kz + <p) \ n _s 

\i(ksin(kz + (i>)-mcos(kz + (f)))/E J' [U > 

Let us analyze the cap region. For small A one has p = i? t \/2A, ^/9zz = oe/V2A,W v ~ —l,a v = ±3/2 and the 
equations (15) are written as 

-i(V2Ad A + -L=(j — B — l/2))v - eu, 
V2A 

-i{V2Ad A --?=(j-B + l/2))u = sv, (18) 
v 2 A 

with B = — 1 ± 3/2, e — ER t . Notice that (18) are similar to the equations for spherical fullerene [5]. At A the 
general solutions (for both u and v) have a power-like form CiA" + ^A - ", a > 1/2. Since the second term cannot 
be normalized, one has to put Ci = 0. 

In the general case, the system (15) can be solved numerically. To this end, the initial conditions are taken to be 
(17) for z » Rf, and the convergence of the wavefunction at z = —Rf is obtained by a variation of the phase <j> up 
to the moment when C2 turns out to be zero within calculation accuracy. The wavefunction is properly normalized, 
which means C = \J DoSo(E) where DoSq{E) = (dE/dk)^ 1 is the density of states in the infinite tube. The results 
of numerical calculations are presented in figs.l and 2. 





FIG. 2. The density of states (per unit area, in arbitrary units) as a function of the energy in the cap region (bottom), near 
the cap (middle) and far from the cap (top). The case of metallic (left) and semiconducting (right) tubes is presented. We take 
l/a=0.9, and the energy E is measured in the units of hVp/Rt- 



4 



Fig. 2 shows the density of electronic states as a function of the energy in three regions: far from the cap, near the 
cap and on the cap for both metallic and semiconducting cases. As is seen, the peaks appear at the energies higher than 
the threshold energy m, where mR t = 0,±1,±2, .. for metallic and mR t — ±1/3, ±2/3, ±4/3, .. for semiconducting 
case. Notice the appearance of secondary (less pronounced) peaks far from the cap. For metallic tube, the constant 
DoS below the threshold energy is found. 




Fig. 3 shows the DoS as a function of the energy and the coordinate. In the cap region, both the peak values and 
the constant DoS (in metallic nanotubes) decrease markedly. At high energies, the harmonic dependence from the 
coordinate is seen in fig. 3 (left), which is in accordance with (17). 

The most interesting finding is the finite values of the DoS near points E — m where the DoS behaves linearly 
in E — m. Moreover, as is seen from figs. 2, 3, the smoothed VHS peaks are shifted to higher energies. Within our 
model the dispersion relation has the form E = Vk 2 + m 2 and, therefore, DoSq(E) should diverge like {E — m)^ 1 / 2 
when the energy approaches the threshold m. One could expect a similar behavior in numerical calculations but this 
was not observed. It should be noted that the only way to avoid the singularity is to put cf> — ±tt/2 in (17). In this 
case both spinor components turn out to be finite at E w m (both u and v become proportional to k). Thus, the 
numerical results clearly indicate that the phase is fixed 

lim 4>{E) = ±vr/2. (19) 

E— >m+0 

Let us proof this analytically. Approaching to the cap, the tube tends to decrease its radius which results in a 
corresponding increase of the effective mass. Using the substitution ip — /(A)V> where /(A) = exp(— ae~ 2A /2) and 
expanding the equations (15) at large A up to the order of e" 2A one obtains 

-i(d A T Ml + e- 2A /2)) (~) =e (l) ' ( 2 °) 

where e = ER t /a, fi — (M2/3 — j)/a. This system is easily reduced to two decoupled equations of the second order 

d 2 A ^ + (Ae- 2A + k 2 )^ = 0, (21) 

where k 2 = e 2 — fi 2 , A = —(p 2 ± /1) where the sign plus (minus) relates to ^ — u (v), respectively. In view of 
M2 = 0, ±1/3, j = ±1,±2, .. and a ~ 1 the constants A for u and v always take opposite signs. For positive A 
(which, for example, relates to v), the substitution x = ^/A^e~ A in (21) leads to the Bessel equation of index in 

x 2 v" + xv' + (x 2 + k 2 )v = 0, (22) 

and for u there appears the modified Bessel equation (with the argument x = iy // \A~^\e^ A ). The solutions should not 
diverge at z — > 00 (x — > 0), so that they are the Bessel and modified Bessel functions of the first kind, respectively. 
For small k (but fc>l) one has in the leading order 
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v w iC v K\ogx — iC v kz, 

u rts iC u n\ogx — iC u kz. (23) 

Substituting this into (20) one obtains the exact relation iC v = C u or iv = u in (17), which determines the limit of 
the phase at E — > m explicitly as <j) — ir/2. Assuming a smooth dependence of the phase <j) on the energy in (17) we 
obtain both the high peak at E > m and the linear dependence of the DoS near E = m. The amplitude of the peak 
is found to grow with increasing z (cf. fig. 3). 

In conclusion, wc have developed a phenomenological model for the description of electronic states in capped 
nanotubes, both semiconducting and metallic. The influence of pentagons in the cap region was taken into account 
via the fictitious compensating gauge fields. We have suggested a single manifold for the description of both the 
nanotube and the cap region. It should be noted that a similar approach was recently used in [12] for the description 
of both curved and defect-free regions in graphcnc. The asymptotical solutions in the tube region were used to obtain 
the initial conditions for the numerical calculations of the density of states. The numerical calculations show that 
the VHS is smoothed out in DoS. This is explained by the fact that the phase of the wavefunction in the tube region 
is determined by the asymptotic solution near the cap. As it was analytically shown, in the cap region the mass 
factor in (20) markedly increases due to geometry. This results in fixation of the wave-function phase which, in turn, 
provides the smoothing of VHS. Thus, the effect of VHS smoothing is geometrical in its origin within our approach. 

The numerical calculations show the finite peaks at energies greater than the threshold energies. The amplitude of 
the peaks grows with z. At threshold points E — m, the DoS is found to be linear in E — m rather than manifests the 
standard singular behaviour. The obtained results are in qualitative agreement with the experimental measurements 
and tight-binding calculations for the closed nanotube performed in [3]. At the same time, near the cap an increase of 
DoS was observed and explained in [3] by the specific distribution of six pentagons on the non-spherical cap. Notice 
that in our model the defects in the cap region are taken into account by effective continuous fields. For this reason, 
the localized states are smeared out and there are no pronounced peaks. Between the peak energies the DoS on the 
cap was found to be reduced in comparison with the tube region [3] . In our model the decrease of the DoS on the cap 
is found at all energies. 
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